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The hadronic hght-by-hght scattering contribution to the muon anomalous magnetic moment can 
be estimated by computing constituent quark loops. Such an estimate is very sensitive to the numer- 
ical values of the constituent quark masses. These can be fixed by computing the hadronic vacuum 
polarization contribution to the muon magnetic anomaly within the same model. In this Letter, we 
demonstrate the stability of this framework against first-order perturbative QCD corrections. 



The measurement of the muon anomalous magnetic 
moment completed about seven years ago [Ij] resulted in 
one of the most interesting controversies in contemporary 
high-energy physics. When the very precise experimental 
measurement 
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is compared to the highly-developed theoretical value [H 
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a 3.2a discrepancy is observed: 
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The size of the discrepancy is tantalizingly close to that 
predicted by the most popular extensions of the Standard 
Model of particle physics, such as Supersymmetry 

The Standard Model value of the muon anomalous 
magnetic moment is obtained by combining QED, elec- 
troweak and hadronic contributions. While the QED and 
electroweak corrections to the muon anomalous magnetic 
moment are well-established, the hadronic contributions 
are notoriously difficult to handle. Indeed, given the 
small value of the muon mass = 105.7 MeV, these 
contributions are sensitive to non-perturbative hadron 
physics. Their calculation from first principles is nearly 
impossible except, perhaps, by lattice field theory meth- 
ods 4, 5]. 

It is useful to distinguish two types of hadronic con- 
tributions to the muon magnetic anomaly - the hadronic 
vacuum polarization and the hadronic light-by-light scat- 
tering. The calculation of the hadronic vacuum po- 
larization contribution uses dispersion relations to con- 
nect it to the integral of the annihilation cross-section 
a{e~^e~ — hadrons) over center-of-mass energies. This 
cross-section has been measured over a wide range of 



energies and is known precisely. The hadronic vacuum 
polarization contribution to the muon anomalous mag- 
netic moment can therefore be calculated in an essentially 
model- independent way 0| ■ 

On the contrary, the hadronic light-by-light scattering 
contribution to the muon anomalous magnetic moment 
can not be directly related to any experimental data. Its 
calculation is therefore bound to rely on model assump- 
tions. Popular models that address the computation of 
the hadronic light-by-light scattering contributions uti- 
lize the existence of two parameters in QCD - the num- 
ber of colors Nc and the mass of the pion m,^. It is 
well-known that in the limit of large Nc or small niTr, low- 
energy QCD simplifies. In the chiral rriTr limit, the 
dominant contribution to the muon anomalous magnetic 
moment comes from loops of charged pions whose inter- 
action with photons is given by scalar QED, up to power 
corrections. In the large-iVc limit, the dominant con- 
tribution comes from approximating the photon-photon 
scattering by single-meson exchanges. 

It is a feature of all existing model calculations of the 
hadronic light-by-light scattering contribution [8.-,12j that 
the chiral enhancement is not effective. A plausible ex- 
planation of this feature is given in Ref. 113], where the 
anatomy of the charged pion loop contribution, includ- 
ing the power-suppressed terms, is studied. The absence 
of a meaningful chiral enhancement leaves us with the 
number of colors Nc as the only parameter to employ. 
Interestingly, the dominance of this parameter allows us 
to extrapolate the calculation from the low-momentum 
region to high-momentum region since the quark contri- 
bution to hadronic light-by-light scattering is enhanced 
in the large- iVc limit. As pointed out in Ref. [12], this 
feature leads to a useful QCD constraint on the hadronic 
light-by-light scattering amplitude. 

It is possible to take these ideas to the extreme and 
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estimate the full hadronic light-by-light scattering con- 
tribution to the muon anomalous magnetic moment from 
constituent quark diagrams alone. This approach has a 
long history. Its early stages are described in Ref. [Tsj . 
The quark loop computation was used to question the 
sign of a^^^ in Refs. [14| . | l 5l| and, gradually, it started 
becoming presented [la, Il6l| as a respectable way to esti- 
mate the hadronic light-by-light scattering contribution 
to the muon magnetic anomaly. 

There are two immediate problems that a computa- 
tion of the hadronic light-by-light scattering through the 
constituent quark loop must address. The first prob- 
lem is that a^^^ depends strongly on the constituent 
quark masses. However, those masses can be estimated 
by requiring that the same model correctly predicts the 
hadronic vacuum polarization contribution to the muon 
anomalous magnetic moment. Since the latter can be 
also obtained in a model-independent way by integrat- 
ing the annihilation cross-section cr(e^e~ — ^ hadrons), 
we get a very useful constraint. Such a procedure was 
suggested in Ref. fl§\ and later adopted in Ref. ^16^. In 
both of these references, subtle details including SU (3)- 
flavor mass differences of constituent quarks were in- 
cluded. It is interesting that this approximate framework 
appears to be remarkably consistent with other com- 
putations based on different models for photon-hadron 
interactions. For example, in Ref. [l6l | the hadronic 
light-by-light scattering contribution was estimated to 
be a'^^^^ = (137^15) ^ 10"^\ in nearly perfect agree- 
ment with the result of Ref. [l^- A more recent esti- 
mate of the hadronic light-by-light scattering contribu- 
tion aji'^i = (106 ± 25) X 10"" was given in Ref. [13. 
It attempts to accommodate many existing estimates of 
this quantity in a self-consistent manner, and is consis- 
tent with the results of Refs. [Tsl. [l6j. 

The second problem of the constituent quark model is 
that it does not possess the correct chiral limit of QCD 
unless quarks directly couple to the Goldstone bosons, 
i. e. the pions. Such extensions of the constitue nt q uark 
model are known and are well-documented |2(]I422 |. In 
the context of the hadronic-light-by-light scattering con- 
tribution to the muon magnetic anomaly, the tt" couples 
to photons through a constituent quark loop. This in- 
troduces a double logarithmic contribution to the muon 
magnetic anomaly, a'^^^ ^ In^ Mq/rriT^ [l^l- This chi- 
ral logarithm cannot be reproduced by a loop of con- 
stituent quarks that couple to photons only. However, 
empirical evidence suggests that the leading logarithmic 
contribution is not numerically important; for example 
see Ref. 24]. Moreover, the logarithmic enhancement be- 
comes ineffective if the numerical values of quark masses 
become comparable to the physical pion mass. As we 
will see in what follows, the quark masses obtained from 
fits to hadronic vacuum polarization are relatively small. 
This provides justification for neglecting diagrams with 
the ttq exchanges. The largest contribution is then due 
to the constituent quark loop. We henceforth focus on 
this term. 



Before discussing the constituent quark loop in de- 
tail, we point out that the good agreement between re- 
sults for aj^j"^' obtained by different groups over many 
years was recently challenged by the results presented in 
Ref. [18,]. In that reference, a computation based on the 
Dyson- Schwinger equation in QCD lead to a result for 
the hadronic light-by-light scattering larger than other 
computations by about a factor of two. The authors of 
Ref. [ill attributed their result to large modifications of 
the quark-photon vertex, when this vertex is computed 
dynamically in nearly full QCD. 

It is important to understand and possibly cross-check 
the results reported in Ref. [ll]. Regardless of the out- 
come, a resolution of this problem will increase our con- 
fidence in the computation of the hadronic light-by-light 
scattering contribution to a^. The need for that is quite 
urgent, given plans to continue the g — 2 experiment at 
Fermi National Accelerator Laboratory, and to reduce 
the experimental error in a^'^P down to 20 x 10"^-'^ [l9| . 
Such an experimental precision approaches the uncer- 
tainty in many existing estimates of the hadronic light- 
by-light scattering contribution. It is definitely much 
smaller than the difference between central values in aj^"^' 
quoted in Ref. and Ref. [Ts^. 

In this Letter, we take an initial step towards a better 
understanding of the hadronic light-by-light scattering 
contribution to the muon magnetic anomaly by calcu- 
lating the QCD radiative correction to the constituent 
quark loop. Since the momentum scales in this problem 
are so low, one may wonder whether or not computation 
of QCD corrections is at all meaningful. The point of 
view that we will adopt here is that a constituent quark 
model with dynamical gluons and pions is a consistent 
field theoretical framework to describe QCD phenom- 
ena at moderate energies [1^. As we explained above, 
there are good reasons to expect the constituent quark 
loop to dominate. Studying constituent quark loop dia- 
grams and QCD corrections to it is also interesting for 
the following reason. Unlike Refs. [l^, [l^ this model 
has two parameters - the constituent quarks mass and 
the value of the strong coupling constant at low energies. 
This allows us to test if QCD dynamics affects ajj'*^' and 
aJj^P in a similar or totally different way. In particular, 
QCD radiative corrections directly check the suggestion 
of Ref. [ll] that large modifications of the photon-quark 
vertex may occur in the hadronic light-by-light scattering 
contribution, but can remain undetected in the hadronic 
vacuum polarization. When calculations of the hadronic 
vacuum polarization and the hadronic light-by-light are 
truncated at the same order in the 0{as) expansion, the 
model gives a unique predictive relation between the two 
quantities. 

Given the approximate nature of our computation, we 
assume that constituent quark masses are significantly 
larger than the mass of the muon and we work to leading 
order in m^/Mq [sl]- We completely neglect the SU{'i) 
mass splittings and we use M„ = Md = Ms = Mq. In 
such an approximation, the hadronic vacuum polariza- 
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tion contribution and the hadronic hght-by-hght scatter- 
ing contribution read [1^ [1^ 



-.hvp _ 



,hlbl 



(- 



2 NrmKQl) 



a 



45Af2 



:C(3) - 



(4) 



In Eq. (gl), (Q^) with n = 2,4 denote sums over the 
electric charges of u, d and s quarks in the appropriate 
powers. Taking the ratio of two contributions in Eq. (j4|), 
we find 
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We use a^JP = 6900 x 10"" and find a'i™ = 148 x 10"", 
in reasonable agreement with previous estimates jl7| . 

We point out that these numerical results for the 
hadronic vacuum polarization and the hadronic light-by- 
light scattering contributions require small quark mass 
values, Mq k, 200 MeV. This makes the physical inter- 
pretation of the parameter Mq somewhat obscure. In- 
deed, it is well understood that the p-meson with the 
mass Mp ~ 770 MeV gives the largest contribution to 



iI^^P. Since the mass scale 2M„ 
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icantly smaller than Mp, the fact that one can fit the 
hadronic vacuum polarization contribution to with 
these mass parameters emphasizes how unphysical they 
are. In fact, the contribution of the p-meson to aJj'^P is 
determined by two independent parameters - the mass 
of the p meson and the p — 7 mixing parameter gp ^ 5. 
However, in the constituent quark approximation, the 
quark mass parameter Mq fits simultanesouly their com- 
bination. One finds Mq ^ Mp/ gp which justifies the low 
value of the quark mass. An apparent relevance of this 
small mass scale for other observables, beyond hadronic 
vacuum polarization, suggests an existence of an infrared 
(below the p mass) duality between hadronic and quark 
contributions, whose origin and precise meaning remain 
unclear at the moment '27]. 

We would like to check the sensitivity of Eq.® to 
QCD radiative corrections. The 0{as) correction to the 
hadronic vacuum polarization contribution can be read 
of from Ref. We find 



it is important to check how such modifications of model 
parameters change the hadronic light- by- light scattering 
contribution. It is clear that if the hadronic light- by- light 
scattering contribution receives corrections that are sim- 
ilar to Eq. (|6]), the prediction shown in Eq. ([5]) is hardly 
affected. 

To check if this is indeed the case, we must compute 
the QCD correction to ajj"^' in Eq.Q. Such a calcula- 
tion involves sixty four-loop diagrams. Representative 
examples are shown in Fig. [1] In general, the computa- 
tion of four-loop diagrams is very difficult, if not impossi- 
ble. However, in our case it is relatively straightforward 
for two reasons. First, because we are interested in the 
computation of the anomalous magnetic moment, no mo- 
mentum is transfered from the magnetic field line to the 
muon line, so that relevant diagrams can be reduced to 
self-energy diagrams. Second, since we are interested in 
the limit Mq 3> trip, we can expand those diagrams in a 
Taylor series in nUp and the incoming momentum of the 
muon line, reducing the self-energy diagrams to four-loop 
vacuum bubbles. Upon expanding, all loop propagators 
have momenta that are parametrically proportional to 
heavy quark masses. This momentum region gives the 
only contribution to the muon magnetic anomaly that 
scales as m^/M^. 

Four-loop vacuum bubble diagrams with a single mass 
MeV is signif- parameter Mq are well studied in the literature. They 
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where Cp = 4/3 is a Casimir invariant of the SU{3) 
color group. It is clear from Eq. ([6]) that the QCD cor- 
rection is significant. For numerical estimates, we take 
as = 0.35, which corresponds to renormalization scales 
of about 1 GeV. We find that the QCD corrections to 
^hvp.NLO amount to sixty percent. While this is a large 
correction, it is not entirely meaningful since we can re- 
define the numerical value of the quark mass Mq to ab- 
sorb it, thereby avoiding any change in aJj^P. However, 



can be evaluated u sing various implementations of the 
Laporta algorithm [29j for solving integration-by-parts 
identities [30|. The master four- loop vacuum bubble in- 
tegrals that we require are also known and can be found 
in Ref. 31]. Performing all the necessary steps, we arrive 
at the following result 
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where Aq = -19/16 + 3/2C(3) and 
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In Eq. ^ a4_5 — Li4^5(l/2) are the values of polylog- 
arithmic functions and ^3^5 are the values of Riemann 
zeta-function at the respective argument. Numerically, 
Ci?Ai/Ao evaluates to 3.851. 

We take the ratio of Eqs. (j6l7p . equate it to the ratio 
of physical quantities and obtain 
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where i?NLO ^ ^hibi.NLO /ghvp.NLO^ rpj^g Qf p^LQ 

a /3. 19\ 45(Q4) 



contributions reads 
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FIG. 1: Examples of diagrams that contribute to higher-order 
corrections to the hadronic light-by-light scattering contribu- 
tion to the muon magnetic anomaly in the constituent quark 
model. 

The effects of corrections are absorbed into the func- 
tion fioLs), which reflects the relative sensitivity of the 
hadronic vacuum polarization and hadronic light-by-light 
contributions to radiative corrections. We find 
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1-1-3.851 X 
l + bmix' 
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This function depends only weakly on the strong cou- 
pling constant, changing from /(O) = 1, to /(I) = 0.8, 
to f{oo) = 0.76. Taking aJ^^P = 6900 x 10"" and 
in the interval [0, ..,7r], we find a|^"^' to be in the range 
= (118 — 148) X 10^^^, consistent with earlier model 
computations ^17]. The QCD effects tend to slightly 
decrease the hadronic light-by-light scattering contribu- 
tions to the muon magnetic anomaly relative to hadronic 
light-by-light scattering, since f{as) < 1. 

In summary, we studied the prediction of the con- 



stituent quark model for the hadronic light-by-light scat- 
tering contribution to the muon magnetic anomaly, in- 
cluding QCD radiative effects. While it is known that 
pious need to be explicitly introduced into the model to 
correctly describe the chiral limit of QCD, we argued 
that, for the muon magnetic anomaly, this issue is nu- 
merically not important. QCD radiative corrections to 
constituent quark loop diagrams provide a tool to diag- 
nose the sensitivity of the vacuum polarization and the 
light-by-light scattering contributions to various aspects 
of QCD dynamics. We demonstrated that the ratio of 
hadronic vacuum polarization and hadronic light-by-light 
scattering contributions to the muon g — 2 is remark- 
ably stable against QCD radiative corrections. Within 
our approach, we are not able to detect a large renor- 
malization of the quark-photon vertex - particular to 
hadronic light-by-light scattering contribution - that is 
claimed to be responsible for the large enhancement of 
^hibi observed in Ref. 18]. We believe that our result 
indicates that the constituent quark model computation 
of the hadronic light-by-light contribution to the muon 
magnetic anomaly is robust. 
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